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DETECTION'  Or  SURFACE  WAVES  FROM  SMALL  EVENTS  AT  TEDESEISMIC  DISTANCES 


ABSTRACT 


A  matched  filter  approach  for  distinguishing  weak  teleseismic 
surface  wave  signals  from  background  noise  is  presented.  The  method 
diacriminatea  against  events  not  located  in  a  particular  source  re¬ 
gion  of  interest  and  provides  estimates  of  magnitude  and  radiation 
pattern,  when  a  number  of  recording  stations  are  available.  Test 
cases  and  typical  results  for  different  source  regions  are  discussed. 

INTRODUCTION 

One  of  the  important  measures  of  source  mechanism  is  the  ex¬ 
citation  of  surface  waves.  This  encompasses  such  estimates  as  total 
energy,  radiation  pattern,  frequency  spectrum,  and  radiation  of  Love 
wave  energy  relative  to  Rayleigh  wave  energy.  At  one  time  or  another 
practically  all  of  these  estimates  have  been  applied  in  the  study  of 
large  magnitude  events.  However,  because  the  surface  waves  from 
small  events  (M  <_  5)  usually  arrive  at  teleaeiamic  distances  with 
amplitudes  at  or  below  the  noise  level,  special  techniques  must  be 
employed  to  extract  the  desired  information.  We  present  here  a 
matched  filter  approach  foi  detecting  Bmall  amplitude  surface  waves, 
estimating  their  total  energy  content,  and  determining  their  radia¬ 
tion  pattern.  Results  for  several  test  cases  and  several  events  are 
shown,  to  demonstrate  the  usefulness  of  this  approach  as  well  as 
some  of  its  limitations. 

METHOD 

Basically  the  matched  filter  approach  amounts  sitqply  to 
searching  a  record  x(t)  for  a  known  waveform  y(t).  In  particular. 


-1  - 


it  is  assumed  that  x(t)  =  ay(t)  +  n(t)  where  a  is  a  constant  and 
n(t)  is  a  random  noise  process.  If  no  further  assumptions  are 
made  with  regard  to  the  nature  of  the  noise  process,  one  can 
determine  the  least-squares  estimate  of  a^  which  minimises  J(a,T  )  = 

—  2 

y  Tx(t  +  t)  -  a  y(t)  ]  where  the  summation  is  over  the  length  of 

y(t).  In  the  test  series  x,  which  is  of  longer  duration  than  y, 

the  lag  T  indicates  at  what  point  in  x  the  comparison  is  begun. 

The  value  of  obtained  by  setting  d  j  =  0  is 

a  a 

a  =  Z  x<t  +  t)  y ( t ) /  Z  y2lt) 

11  (1) 

Thus,  the  matched  filter  in  this  case  is  simply  the  waveform  y(t) 
and  the  matched  filter  output  at  lag  t  is 

Z  x (t  +  t)  y(t)  =  Cxy  (t)  (2) 

The  coherency  at  lag  i  ir.  given  by 

C(t)  =  Z  x(t  +  t)  y  ( t )  /  [  Z  x2(t  +  t)  •  Z  y2(t)]1? 

t  t  t  (3) 

It  always  is  bounded  by  -1  £  C(t  )  £  1.  The  maximum  in  the  envelope 

C(t  )  occurs  at  the  value  of  t  where  x  and  y  match  best  in  the 

least-squares  sense.  This  C  (  t  )  is  the  correlation  coefficient. 
^  max 

In  the  present  application  we  have  used  this  simple  least 
squares  approach.  Estimates  of  ja,  C  ( T  ) ,  and  C(T  )  calculated 
in  this  way  allow  us  to  study  test  cases  and  unknown  events  with 
a  minimum  of  assumptions. 
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A  further  refinement  in  the  method  can  be  made,  however,  if 
-he  noiae  is  a  stationary  random  process  whose  correlation  matrix 
R  can  be  calculated.  If  the  sampled  x(t)  and  y(t)  are  considered 
as  row  vectors,  x  and  y,  then  an  estimate  of  £  is  given  by 

a  -  xR_1  y1/  yR'1  y‘  (4) 

where  the  prime  denotes  transpose.  This  is  the  maximum  likelihood 
estimate  of  h.  In  this  case  R  1  y*  is  the  filter.  (Details  of  the 
statistical  analysis  leading  to  these  estimates  are  given  in  Chapter 
II  of  Reference  1) . 

This  modified  approach  may  allow  us  to  do  a  better  job  in 
detecting  events,  since  it  accounts  for  the  noiae  as  well  as  tne 
known  waveform.  However,  we  have  not  yet  investigated  its  use¬ 
fulness  sufficiently  to  determine  whether  the  extra  computation 
required  to  calculate  the  correlation  matrix  of  the  noise  is  war¬ 
ranted. 

In  the  least-squares  approach  we  chose  the  following  decision 
criteria  to  determine  whether  or  not  a  signal  was  present! 

1.  A  relative  maximum  in  the  envelope  of  C(t  )  must  fall  in 
the  expected  time  window  based  on  travel-time  information. 

2.  This  value  of  C  (t  )  must  be  greater  than  or  at  least 

max 

comparable  to  typical  peak  vilues  in  the  envelope  of  C  (T)  outside 
the  expected  time  window. 

3.  An  arrival  must  be  detected  within  the  proper  expected 
time  window  at  each  of  several  stations. 

If  all  three  of  these  criteria  were  met,  then  the  signal  was 
considered  to  be  present. 


IMPLEMENTATION 


The  matched  filter  approach  just  discussed  requires  that  one 
knows  the  desired  or  expected  waveform  in  order  to  search  effect¬ 
ively  for  that  waveform  in  a  noisy  record.  in  order  to  obtain  a 
suitable  "expected"  waveform  and  at  the  same  time  assure  that 
propagation  effects  would  be  properly  accounted  for,  we  chose  the 
surface  wave  from  a  larger  event  in  the  source  region  of  interest 
as  the  filter  y(t).  Thus,  no  matter  how  complicated  the  paths 
of  propagation  are,  the  paths  for  other  events  in  the  same  source 
region  will  be  nearly  coincident  with  those  of  the  larger  event 
which  produced  y(t) ,  so  the  only  major  differences  in  recorded 
waveform  will  be  source  differences.  This  situation  is  illustrated 
schematically  in  Figure  1. ,  which  shows  the  surface  wave  trajectories 
appropriate  tor  a  source  region  in  the  ocean  and  a  continental  re¬ 
cording  station. 

This  coincidence  of  travel  path  simply  assures  that  the 
earth's  transfer  function  is  the  same  for  all  events  occurring 
in  that  source  region;  it  does  not  mean  that  dispersion  during 
propagation  to  teleseismic  distances  is  unimportant.  To  the  con¬ 
trary,  it  is  this  physical  effect  (diepersion)  which  is  vital  to 
the  success  of  the  technique,  since  it  transforms  the  source  pulse 
into  an  oscillatory  signal  which  is  of  long  dura.tion.  As  the 
duration  of  the  signal  increases,  random  correlations  of  the  signal 
with  the  noise  become  poorer,  with  the  result  that  the  "false  alarm" 
level  is  reduced.  From  this  standpoint  the  larger  the  epicentral 
distance  the  better  the  method  should  work,  since  the  signal  in¬ 
creases  in  duration  the  larger  the  distance.  However,  the  energy 
density  of  the  signal  is  reduced  due  to  dispersion  and  attenuation 
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during  propagation,  so  that  the  signal  to  noise  ratio  is  reduced 
with  increasing  distance.  Therefore,  for  a  given  transfer  function 
and  noise  level,  one  would  expect  there  to  be  an  optimum  observing 
distance  for  deteeMna  events  using  the  matched  filter.  But,  be¬ 
cause  of  the  different  nomc  among  available  stations  and 

the  variety  of  propagation  paths  to  these  sratrons,  it  does  not 
appear  to  be  practical  to  determine  these  optimum  observin'  dis¬ 
tances  experimentally,  what  we  can  estimate,  however,  is  the 
minimum  S/N  at  which  the  technique  will  detect  observed  or  expected 
signal  waveforms. 

Because  of  dispersion  and  frequency-dependent  attenuation  with 
distance,  the  signal  waveform  will  in  general  be  different  at  each 
station  and  thus  a  different  filter  for  each  station  and  source 
region  must  be  used.  However,  this  presents  no  problem  since  a 
single  large  event  in  a  given  source  region  provides  for  each 
station  a  filter  suitable  for  detecting  other  events  located  in 
that  region,  as  well  as  the  temporal  pattern  of  arrival  times  at 
the  station*.  This  pattern  and  the  fact  that  the  filters  are  dif¬ 
ferent  for  each  station  can  be  used  to  advantage  in  discriminating 
against  events  outside  the  source  region  of  interest,  since  one  or 
more  of  the  decision  criteria  listed  earlier  will  not  be  satisfied; 
in  particular,  the  correlation  peak  values  will  be  degraded,  and 
the  arrival  times  of  these  peaks  will  not  produce  the  appropriate 
pattern  of  arrival  times  at  the  observing  stations. 

When  a  small  event  has  been  detected  by  this  method,  then 
the  question  arises  as  to  what  details  we  can  extract  about  the 
source.  The  estimates  of  A  provide  a  convenient  means  of  comparing 
the  magnitude  of  the  small  event  to  that  of  the  larger  one,  since 
is  an  estimate  of  the  ratio  of  the  energy  in  the  small  event  to 
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that  in  y(t) .  Some  care  must  be  taken,  however,  if  the  reference 

event  exhibits  a  strong  radiation  pattern  which  is  different  from 

that  of  the  smaller  event.  In  fact,  the  radiation  pattern  of  the 

■mall  event  relative  to  the  reference  event  can  be  obtained  by 

plotting  a/C  y^ ( t)  vs.  source-station  azimuth, 
t 

A 

Comparisons  of  a  vs.  distance  along  a  given  azimuth  allow 
us  to  estimate  roughly  the  excitation  spectrum  of  the  small  event 
compared  to  the  larger  one,  since  the  attenuation  with  distance  is 

A 

frequency-dependent,  which  implies  that  a  Bhould  vary  with  distance 
in  a  way  which  reflects  the  shape  of  the  spectrum  cf  the  small  event 
compared  with  that  of  the  reference  event.  Only  if  their  epectra  are 
identical  in  shape  will  a  be  invariant  with  distance.  Since  this  at¬ 
tenuation  factor  is  of  the  form  exp  [-g(u>)&]  ,  where  g(m  )  is  an  in¬ 
creasing  function  of  frequency  and  L  is  epicentral  distance,  a  will 
decrease  with  distance  if  the  spectrum  of  the  small  event  ia  peaked 
at  a  higher  frequency  than  the  reference  event.  A  method  for  cal¬ 
culating  the  spectrum  of  the  small  event  relative  to  that  of  the 
reference  event  is  given  in  Appendix  I. 

Thus  potentially  the  matched  filter  method  permits  us  to: 

1.  detect  the  surface  waves  from  a  source  region  of  interest, 
while  rejecting  events  outside  the  source  region? 

2.  make  magnitude  estimates  for  small  events  relative  to 
the  reference  event? 

3.  outline  radiation  patterns  for  small  events? 

4.  estimate  spectral  shape  relative  to  the  reference  event. 

Computer  programs  were  written  to  implement  the  method  along 
the  lines  discussed  above. 
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TEST  CASES 


In  order  to  teat  the  effectiveness  and  sensitivity  of  the 
matched  filter  approach  for  surface  waves,  several  test  cases 
were  investigated. 

Case  I  -  Known  signal  added  to  actual  noise  at  different  signal -to- 
noise  (S/N)  ratios 

Figure  2  shows  the  results  of  adding  the  signal  shown  at  the 

top  (an  actual  oceanic  Rayleigh  wave  train  recorded  at  4686  km 

distance)  to  the  noise  trace  at  a  number  of  arrival  times,  with  the 

indicated  signal/noise,  S/N,  and  using  the  signal  as  the  matched 

filter.  We  define  S/N  as  m|x  |  y  ( t]j  /RMS  [n(t)]  .  The  number  in 

parentheses  is  an  alternate  value  of  signal/noise  (S/N)  given  by 

RMS  [y(t) ]  /RMS  fn(t) ]  .  The  values  of  j  •  S  are  printed  next  to 

C  trace,  and  the  correlation  coefficients  CC  next  to  the  C  trace, 
xy 

The  arrows  indicate  the  beginning  time  for  the  signal  in  each  in¬ 
stance.  In  this  example  the  signal  is  detected  for  s/n  (S/n)  as 
low  as  .35  (.15),  where  visual  detection  is  impossible.  Notice 
also  that  the  a's  are  approximately  in  proportion  to  the  s/N  used 
as  expected. 

The  reason  the  results  appear  to  be  as  good  or  better  for 
S/N  *  .35  compared  to  S/N  =  .5  is  that  the  noise  level  was  not 
uniform  over  the  entire  noise  trace  used  to  establish  the  RMS  of 
the  noise  and  the  .35  signal  was  buried  at  a  position  where  the 
actual  noise  level  over  the  filter  length  was  low  relative  to  that 
where  the  .5  signal  was  buried. 


Figure  3  shows  the  results  when  a  different  signal  (  a 
synthetic  seismogram  calculated  for  an  assumed  dispersion  curve) 
is  added  to  the  noise  at  different  levels.  The  format  is  the 
same  as  in  Figure  2.  This  signal  waveform  approximates  that  for 
a  purely  continental  propagation  path  with  the  dispersion  shown 
An  Figure  4.  We  can  detect  this  signal  down  to  a  S/N  of  about 
.4  (  . J  4) ,  which  is  again  considerably  below  the  threshold  of 
visual  detection.  This  noise  sample  was  identical  to  that  used 
in  Case  1,  so  the  same  comments  hold  with  regard  to  how  well  the 
signal  is  detected  at  different  positions.  These  cases,  along  with 
results  for  other  test  cases,  suggest  that  if  the  signal  buried  in 
the  noise  is  identical  in  waveform  to  the  reference  event,  it  can  be 
detected  when  the  S/N  is  above  about  .35  (.15)  and  that  the  a's  pro¬ 
vide  reliable  measures  of  the  relative  amplitudes  of  the  two  events 
down  to  a  s/N  of  about  .5.,  It  should  be  noted  that  the  indicated 
a  values  include  a  gain  factor  introduced  in  constructing  a  partic¬ 
ular  S/N  level  in  the  x(t)  trace  using  y(t)  and  a  given  n(t). 

Case  III  -  Effects  of  varying  eplcentral  distance 

In  order  to  test  the  sensitiveness  of  the  method  to  pertur¬ 
bations  in  source  distance,  we  wrote  a  computer  program  to  syn¬ 
thesize  the  surface  train  to  be  expected  at  any  distance  for  a 
given  phase  velocity  dispersion  curve  and  amplitude  spectrum 
(see  Appendix  II) .  Using  the  dispersion  curve  shown  in  Figure  4 
and  the  amplitude  spectrum  labeled  Aq  in  Figure  7,  we  synthesized 
seismograms  for  distances  in  the  range  from  3000  to  4000  km,  in 
50-km  increments.  Figure  5  shows  some  of  these  synthetic  records 
and  the  corresponding  matched  filter  results,  when  the  seismogram 
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•t  3000  km  was  used  as  the  reference  in  each  case  and  no  noise 
was  present.  The  values  of  a  and  CC  are  a  measure  of  the  effect 
of  distance  perturbations  on  tho  matched  filter  response.  Sur¬ 
prisingly  perhaps,  the  method  appears  to  be  rather  insensitive  to 
changes  in  source  distance;  this  is  indicated  in  Figure  6  where 
the  correlation  coefficients  are  plotted  vs.  differences  in  epi- 
central  distance,  since  the  group  velocity  window  and  amplitude 
spectra  were  kept  fixed  in  this  case  the  values  of  a  and  CC  are 
identical,  lhere  is  the  expected  decrease,  but  the  correlation 
coefficient  is  above  .8  for  distance  differences  as  large  as 
1000  km.  It  should  be  pointed  out  that  this  result  is  not  inde¬ 
pendent  of  the  assumed  dispersion  and  amplitude  spectrum.  That  is, 
a  given  perturbation  in  distance  will  produce  a  change  in  signal 
waveform  which  depends  on  both  the  dispersion  in  the  source  region 
and  the  amplitude  spectrum.  Nevertheless,  this  case  is  probably 
sufficiently  representative  to  infer  that  the  method  does  not  de¬ 
pend  critically  on  the  closeness  cf  the  reference  source  to  the 
source  of  interest  insofar  as  the  correlations  are  concerned.  This 
was  borne  out  in  the  analysis  of  actual  events.  However,  the  rela¬ 
tive  travel-time  pattern  would  atill  require  the  events  to  be  rather 
close  together,  since,  for  example,  a  100-km  difference  could  result 
in  relative  arrival  time  shifts  among  stations  of  as  much  as  25-30 
seconds. 

Case  IV  -  Effects  of  amplitude  spectra  on  matched  filter  response 
Differences  in  the  amplitude  spectrum  of  the  reference  event 
compared  to  the  test  event  can  be  expected,  since  in  general  the 
lower  the  magnitude  the  higher  the  frequency  of  the  peak  in  the 


amplitude  spectrum  observed  at  a  given  distance.  Thus,  when  we 
look  for  small  events  using  large  ones  we  almost  certainly  are 
not  dealing  with  spectra  which  are  aimilar  either  in  shape  or 
position  in  frequency.  To  test  these  effects  we  computed  the 
matched  filter  response  for  the  synthetic  seismograms  corresponding 
to  the  anplitude  spectra  shown  in  Figure  7.  The  dispersion  used 
was  the  same  as  for  Case  III;  the  distance  was  taken  as  3000  km, 
and  the  reference  spectrum  is  the  one  labeled  AQ  in  Figure  7.  The 
matched  filter  results  are  shown  in  Figure  8.  In  this  case  the 
effect  on  both  a  and  the  correlation  coefficient  is  more  pronounced, 
as  shown  in  Figure  9.  One  reason  that  a  and  CC  do  not  coincide  is 
that  only  energy  in  a  fixed  group  velocity  window  is  included  in  the 
synthesis,  so  the  relative  energies  may  differ  from  one  calculation 
to  another.  In  any  case  our  estimates  of  &  and  CC  are  apt  to  be 
lowered  if  the  excitation  spectra  are  different.  On  this  basis  it 
would  be  desirable  to  use  smaller  magnitude  eventB  for  reference 
sources,  but  if  this  is  done  then  the  reference  signal  itself  is 
degraded  by  noise,  we  conclude  from  visual  examination  of  the 
teleseismic  Rayleigh  waves  from  a  number  of  events  that  earthquakes 
in  the  body  wave  magnitude  range  5. 5-6.0  usually  have  a  high  enough 
s/n  to  be  usable  as  reference  events.  Some  lower  magnitude  earth¬ 
quakes  also  exhibit  large  enough  surface  waves  to  be  used  for 
reference. 
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results 


In  this  section  we  present  sose  of  the  results  of  applying 
the  matched  filter  approach  discussed  above  to  actual  data.  These 
results  are  intended  primarily  to  demonstrate  the  practical  appli¬ 
cability  of  the  method,  and  no  claims  are  made  about  a  general 
minimum  threshold  of  detection,  although  events  as  small  as  mag¬ 
nitude  4.1  were  observed  at  taleseismic  distances  in  some  instances* 
Many  more  small  events  must  be  analysed  before  the  performance  limits 
of  the  matched  filter  approach  can  be  adequately  assessed. 

Results  for  three  different  source  regions  -  Aleutians  (Am- 
chitka) .  Kamchatka,  and  Hawaii  -  are  presented.  The  observation 
stations  are  LRSM  stations  or  VELA  Observatories  in  North  America. 

The  epicenter  information  for  the  three  sets  of  events  we  discuss 
is  given  in  Table  I. 

The  organisation  in  the  figures  to  follow  conforms  to  the 

standard  format  adopted  for  displaying  all  matched  filter  results 

for  actual  events.  That  is,  a  station  set  consists  of  four  traces; 

the  top  trace  is  the  signal  y(t),  the  second  trace  is  the  observed 

seismogram  x(c),  the  third  trace  is  C  (t  ).  end  the  bottom  trace  is 

a 

C(  t).  The  legend  at  the  left  gives  values  of  distance,  a,  and  CC 
as  well  as  station  and  event  identification.  The  arrows  indicate 
the  location  of  the  expected  signal  arrival  at  each  station.  This 
ie  where  the  peak  in  the  matched  filter  output  should  occur;  that  is, 
the  beginning  point  of  y(t)  in  x(t).  The  interval  between  each  time 
mark  is  50  seconds  unless  otherwise  indicated, 

ALEUTIAN  SOURCE  REGION 

Comparison  of  an  explosion  and  an  earthquake  of  comparable  macmitude. 

The  explosion  was  L0NGSH0T  (body  wave  magnitude  5.97);  the 
earthquake  was  a  magnitude  5.9  event  located  within  a  degree  of  the 


11 


LONG SHOT  sit*  (•••  Table  I) .  At  each  station  ths  surface  waves 
from  the  22  November  1965  earthquake  were  chosen  as  the  reference 
signal,  and  the  LONG SHOT  seismograms  were  scanned  for  surface  waves 
with  the  earthquake  signal  used  as  the  matched  filter,  some  of  the 
results  for  Rayleigh  waves  are  shown  in  Figure  10. 

It  is  seen  that  the  Rayleigh  wave  was  clearly  detected  at 
all  but  one  station.  The  Rayleigh  wave  amplitude  of  LONGSHOT  rel¬ 
ative  to  that  of  the  magnitude  5.9  earthquake,  however,  is  below  .12 
everywhere  in  the  distance  range  analyzed  (3500-7000  km).  This 
means  that  the  surface  wave  magnitude  for  LONGSHOT  is  only  about  5.0 
if  the  surface  wave  magnitude  of  the  reference  earthquake  is  equiva¬ 
lent  to  its  body  wave  magnitude.  It  also  means  that  the  r/p  exci¬ 
tation  for  the  explosion  is  nearly  an  order  of  magnitude  smaller 
than  the  R/P  excitation  for  an  equivalent  size  earthquake. 

The  results  in  Figure  10  are  for  the  vertical  components  only. 
We  did  the  same  analysis  for  the  radial  components  and  got  essentially 
the  same  results.  Therefore,  these  are  not  shown. 

We  also  tried  different  signal  durations  for  y(t)  and  found 
that  the  signals  shown  in  Figure  10  gave  a  better  result  (correlation' 
than  did  longer  signals  which  included  part  of  the  coda  of  the  earth¬ 
quake.  This  probably  reflects  the  situation  that  in  using  the  longer 
filter  more  noise  was  included  in  the  correlation  interval  for  both 
y(t)  and  x(t)  without  a  proportional  increase  in  signal  energy  in¬ 
cluded  in  the  interval. 

Since  the  22  November  event  also  generated  strong  Love  waves, 
we  used  them  as  a  matched  filter  to  search  the  LONGSHOT  seismograms 
for  Love  waves  in  the  same  way  we  searched  for  Rayleigh  waves.  We 
used  the  transverse  components  throughout,  with  the  possible  excep¬ 
tion  of  two  stations  wo  were  not  able  to  detect  any  Love  waves  from 
LONGSHOT.  Therefore,  in  Figure  11  we  show  the  results  for  only 
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three  atationa  including  the  two  where  e  Love  wave  arrival  might 
be  indicated.  The  bottom  set  in  Figuis  11  is  representative  of 
the  atationa  where  nothing  was  detected.  The  indicated  values 

A 

of  a  in  Figure  11  give  an  upper  bound  of  about  .01  for  the  ampli¬ 
tude  of  Love  waves  for  LONG SHOT  relative  t)  the  Love  waves  for  the 
22  November  earthquake.  This  result  is  not  surprising,  since  an 
examination  of  the  short  period  seismograms  for  LONG  SHOT  rev* als 
that  it  was  an  explosion  which  excited  shear  waves  poorly  or  not 
at  all. 

in  Figure  12  is  shown  a  plot  of  &  vs.  distance  for  the 
LONG SHOT  Rayleigh  waves  observed  along  the  azimuth  shown  in  Figure 

A 

16.  The  observed  values  of  a  decrease  with  distance  which  indicates 
that  the  spectrum  of  LONGSHOT  Rayleigh  waves  is  peaked  at  higher 
frequencies  than  is  the  22  November  earthquake  spectrum  (see  Appen¬ 
dix  I).  In  all  cases  the  values  of  &  are  below  .12  indicating,  as 
mentioned  before,  a  smaller  Rayleigh  wave  magnitude  for  LONGSHOT. 

Comparison  of  smaller  magnitude  earthquakes  with  the  22  November 
1965  earthquake. 

In  order  to  concentrate  on  the  smaller  magnitude  events  in 
the  Amchitka  area,  we  analyzed  only  a  few  events  for  intermediate 
magnitudes  of  approximately  5,  ueing  only  a  few  stations  to  verify 
that  the  Rayleigh  waves  could  be  detected.  The  Rayleigh  waves  from 
earthquakes  of  magnitude  greater  than  5  from  this  region  can  be 
clearly  seen  visually  on  the  film  records,  so  the  matched  filter 
analysis  is  not  needed  for  detection  of  these  events.  Figure  13 
shows  the  matched  filter  output  for  a  single  station  for  the  mag¬ 
nitude  4.9  event  listed  in  Table  I.  The  surface  wave  is  clearly 
detected. 


13 


In  Figure  14  we  show  the  natched  filter  results  for  the 
magnitude  4.0  earthquake  (20  February  iy&6)  listed  in  Table  I. 

Notice  that  the  source  depth  for  this  event  was  48  km  which  means 
that  the  surface  wave  excitation  i*>  necessarily  smaller  than  if 
the  same  source  were  located  closer  to  the  surface.  From  Figure  14 
it  appears  that  this  event  was  barely  detected  at  several  stations 
but  not  detected  at  every  station.  The  values  of  &  average  about 
.01  which  is  approximately  the  value  expected  for  a  magnitude  dif¬ 
ference  of  2.  Certainly  this  event  is  at  cr  near  the  lower  limit 
of  detectability  using  the  matched  filter  approach. 

From  the  limited  analysis  we  have  done,  it  would  appear  that 
we  can  detect  earthquake  events  with  magnitudes  approaching  4  from 
the  Amchitka  regicn.  certainly  the  threshold  is  below  magnitude  5. 
However,  more  small  events  must  be  analyzed  before  a  lower  magnitude 
threshold  for  detection  can  be  firmly  established. 

KAMCHATKA  SOURCE  REGION 

We  present  in  this  section  results  for  the  5  events  located 
in  Kamchatka  which  are  listed  in  Table  I.  In  thene  cases  the  ob¬ 
serving  distances  were  in  the  range  from  5300-8000  km  which  overlaps, 
but  extends  further  than,  the  distances  investigated  for  the  Aleutian 
events.  We  consider  only  Rayleigh  waves  and  use  the  Rayleigh  waves 
from  the  29  January  1965,  magnitude  5.8  event  as  the  matched  filter 
in  all  cases.  Note  in  Table  I  that  all  these  events  are  located 
within  2-3  degrees  of  each  other. 

Analysis  of  the  magnitude  5.3  event  (20  April  1965) 

Shown  in  Figure  15  are  some  of  the  matched  filter  results  for 
this  event.  The  event  clearly  was  detected  at  all  but  2  of  the  20 
stations  analyzed.  These  two  HL2ID  and  blwv,  aro  included  in  Figure  15. 
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They  possibly  define  a  node  on  the  radiation  pattern  since  they  are 

aligned  approximately  on  a  single  azimuth  from  the  source.  Mote 

that  at  some  stations  the  correlation  coefficient  for  the  arrival  is 

quite  high  (  >  .6)  .  Figure  16  is  a  map  showing  the  observed  values 
*  -1 

of  (a)  ,  that  is,  the  amplitude  of  the  largn  event  relative  to  the 

small  one  at  the  various  stations,  with  the  exception  of  the  vicin¬ 
ity  of  the  suspected  nodal  line,  the  values  of  (4i)”^  are  consistent 
in  each  geographical  area  where  several  stations  are  clustered.  For 
the  given  USC&GS  body  wave  magnitudes  the  value  of  £  in  this  case 
should  be  about  .32  or  equivalently  its  inverse  should  be  about  3.1, 
a  somewhat  lower  value  than  actually  observed  at  most  stations.  This 
would  suggest  that  the  magnitude  of  the  20  April  1965  event  is  slightly 

below  5.3.  However,  it  must  be  remembered  that  the  higher  values  of 
»  -1 

<a)  observed  could  also  result  from  a  different  radiation  pattern 
or  an  excitation  spectrum  peaked  at  higher  frequencies  than  the  ref¬ 
erence  event.  Therefore,  the  observed  values  are  probably  reasonable 
for  a  magnitude  5.3  event. 

Analysis  of  the  magnitude  5.1  event  (6  July  1965) 

In  Figure  17  are  shown  some  of  the  matched  filter  results 
for  this  event.  It  is  clear  from  the  high  values  of  the  correlation 
coefficient  that  the  event  was  positively  detected.  Moreover,  the 
values  of  a  are  also  large,  and  one  can,  in  fact,  see  the  signal 
above  the  noise  on  most  of  the  seismograms  in  Figure  12.  visuel 
picks  of  the  amplitude  ratio  for  this  event  consistently  gave  larger 
values  of  a  than  those  machine  picked.  Thus,  the  6  July  1965  event 
seems  to  be  .3  or  less  smaller  in  magnitude  than  the  29  January  1965 
event,  which  means  that  the  estimated  surface  wave  magnitude  of  the 
6  July  1965  event  is  greater  than  5.5.  Moreover,  comparison  of 
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Figures  15  and  17  indicates  that  the  6  July  1965  event  was  definitely 
larger  than  the  magnitude  5.3  event  on  20  April  1965.  Therefore,  the 
body  end  surface  wave  magnitude  estimates  for  the  6  July  1965  event 
differ  by  at  least  0.4,  with  the  surface  wave  magnitude  the  larger. 
Analysis  of  the  magnitude  5.0  event  (14  February  1965) 

Shown  in  Figure  18  are  some  of  the  filter  results  for  this 
event.  The  method  apparently  detects  the  event  but  at  a  time  30-50 
seconds  earlier  than  expected  at  each  station.  This  probably  results 
from  the  fact  that  this  event  was  nearly  4  degrees  closer  t^an  the 

A 

reference  event.  The  a's  are  consistan-  <tu6  average  around  77  im¬ 
plying  that  the  surface  wave  magnitude  of  this  event  is  slightly 
more  than  one  unit  below  that  for  the  reference  event,  that  is, 
somewhat  below  4.8.  Since  the  excitation  spectrum  of  the  small  event 
probably  peaks  at  a  nigher  frequency  than  that  of  the  reference  event, 

A 

it  is  expected  (Appendix  I)  that  the  matched  filter  estimate  of  a, 
hence  magnitude,  will  decrease  with  increasing  epicentral  distance. 
Thus,  the  results  for  this  case  are  reasonably  consistent  with  the 
body  wave  magnitude  of  5.0. 

Analysis  of  the  magnitude  4.4  event  (21  January  1965) 

Shown  in  Figure  19  are  some  of  the  matched  filter  results 
for  the  21  January  1965  overt.  It  should  be  noted  that  this  was  a 
deep  focus  event  (119  km) ,  hence  equivalent  to  a  smaller  mangitude 
shallow  event  as  far  as  Rayleigh  wave  excitation  is  concerned.  It 
appears  that  this  event  was  barely  detected  by  the  method  with  es- 

A 

timates  of  a  around  .03.  This  value  gives  a  magnitude  estimate  of 
about  4.2  for  the  21  January  event,  which  seems  large,  considering 
the  source  depth  of  119  km.  In  any  case  this  event  appears  to 
define  the  threshold  of  detectability  for  the  Kamchatka  region. 
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HAWAIIAN  SOURCE  REGION 

Results  for  two  smell  events  «r*  presented  for  this  source 
region.  The  magnitudes  of  these  events  are  4.4  and  4.1  respectively, 
while  the  reference  event  is  the  magnitude  5.3  event  (11  October  1964) 
listed  in  Table  I. 

Analysis  of  the  magnitude  4.4  evc..t  (7  January  1964) 

Shown  in  Figure  20  are  matched  filter  results  for  the  7 
January  1964  event  at  several  stations.  Tne  event  was  apparently  de¬ 
tected,  but  the  amplitude  estimates  are  quite  small  (around  .02)  a 
factor  of  5  smeller  than  expected  if  the  event  is  truly  of  magnitude 
4.4.  The  large  oscillatory  output  later  in  the  records  of  C  and  C 
represents  the  correlation  of  the  reference  event  with  another  large 
event  located  some  7000  km  more  distant  in  New  Britain.  Thia  cor¬ 
relation  is  of  long  duration  because  the  more  distant  event  has  dis¬ 
persed  a  great  deal  more  in  propagating  the  additional  7000  km  and 
the  filter,  therefore,  "sees"  a  signal  for  a  long  time  in  x(t).  Thia 
example  provides  another  illustration  of  the  fact  that,  the  method  is 
weakly  dependent  on  the  separation  of  the  sources  and  that  arrival 
time  patterns  must  be  used  as  additional  constraints  if  we  wish  to 
detect  only  events  from  a  particular  source  region  of  interest. 

Analysis  of  the  magnitude  4.1  event  (13  August  1964) 

In  Figure  21  are  shown  the  resultB  for  the  shallow  focus 
(11  km)  magnitude  4.1  event  of  13  August  1964.  This  event  is  appa¬ 
rently  detected  at  moat  of  the  stations.  However,  the  corresponding 
correlation  coefficients  and  the  estimates  of  h  are  small-  The 
apparent  magnitude  of  the  small  event  is  about  1.6  smal  ’  an  the 

reference  event.  However,  the  filtering  effect  of  oceanic  jpagation 

of  Rayleigh  waves  will  greatly  attenuate  periods  shorter  than  about 
15  seconds  and  as  the  peak  in  excitation  shifts  to  shorter  periods  for 
the  smaller  magnitude  events,  this  effect  should  cause  very  signifi¬ 
cant  reductions  in  surface  wave  energy,  hence  reduce  the  matched  filter 
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estimates  of  surface  wave  magnitude  at  teleseismic  distances. 
Nevertheless,  we  have  been  able  to  “see"  events  of  magnitude  4.1 
and  4.4  for  the  Hawaiian  source  area. 

DISCUSSION 

In  the  foregoing  sections  we  have  demonstrated  experi¬ 
mentally  that  the  surface  waves  from  small  events  can  be  detected 
at  teleseismic  distances  using  a  matched  filter  method  and  that 
reasonable  estimates  of  some  source  characteristics  can  be  made. 

What  has  not  been  demonstrated  is  the  minimum  threshold  of  detection 
for  the  method.  The  underlying  difficulty  which  alwavs  can  affect 
the  estimated  threshold  seriously  is  that  the  body  wave  magnitudes 
used  as  an  absolute  standard  for  comparing  events  may  be  badly  in 
error.  A  more  reliable  absolute  measure  of  source  strength  is  needed 
to  evaluate  the  performance  of  the  matched  filter  in  detecting  actual 
everts. 

One  of  the  most  useful  features  of  the  matched  filter  method  is 
that  it  permits  one  to  make  meaningful  surface  wave  magnitude  estimates 
even  when  the  signal  to  noise  level  is  so  low  that  one  cannot  visually 
detect  the  presence  of  the  signal.  This  means  that  the  matched  filter 
method  can  be  used  to  obtain  reliable  surface  wave  magnitude  estimates 
for  events  so  small  that  meaningful  AR  estimates  of  magnitude  (Ref¬ 
erence  4)  cannot  be  obtained. 

Several  possible  improvements  in  the  matched  filter  approach 
can  be  made.  One  was  discussed  in  the  Method  section,  that  of  in¬ 
cluding  the  noise  correlation  matrix  to  obtain  maximum-likelihood 
estimates  of  A  instead  of  least  squares  estimates.  Another  is  to 
use  empirically  derived  scaling  functions  to  adjust  the  spectrum  of 
the  reference  event  to  match  that  generally  expected  for  a  small 
magnitude  event?  this  should  improve  the  correlation  between  the 
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reference  and  teat  traces.  Another  is  to  phase  equalise,  and  sum 
over  an  array  of  stations  before  applying  the  matched  filter. 

In  fact  there  is  no  reason  why  one  should  not  make  use  of  all  the 
available  noise  suppression  techniques  before  applying  the  matched 
filter;  this  could  in\prove  the  detection  capability  of  the  matched 
filter  significantly. 

The  indirect  method  presented  in  Appendix  I  for  estimating 
the  spectrum  of  the  small  event  has  not  been  investigated  experi¬ 
mentally  in  more  than  the  rudimentary  way  done  for  LONGSHOT.  Since 
being  able  to  determine  the  frequency  spectrum  of  such  small  events 
would  be  of  very  great  value  in  studying  their  source  mechanisms, 
the  practical  applicability  of  the  indirect  method  should  be  ex¬ 
plored  further. 

CONCLUSIONS 

On  the  basis  of  this  investigation  we  conclude  that: 

1.  In  principle  the  matched  filter  approach  can  be  used: 

(a)  to  detect  weak  surface  wave  signals  for  signal  to 
noise  ratios  as  low  as  about  0.35: 

(b)  to  make  reliable  relative  magnitude  estimates  for 
S/N  as  low  as  .5; 

(c)  to  determine  radiation  pattern  relative  to  a  reference 

event; 

[1/  to  estimate  the  general  shape  of  e  small  event's  am¬ 
plitude  spectrum  relative  to  that  of  the  reference  event. 

2.  The  results  presented  for  a  number  of  actual  events  have 
demonstrated  the  practical  applicability  and  usefulness  of  the  matched 
filter  approach  in  studying  weak  teleseismic  surface  waves. 
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3.  The  matched  filter  used  in  conjunction  with  the  surface 
wave  arrival  time  patterns  appropriate  for  different  source  regions 
of  interest  allows  one  to  discriminate  against  all  but  those  events 
in  the  particular  source  region  of  interest. 

4.  Rayleigh  waves  produced  by  the  LONGSHOT  explosion  wure 
detected  at  teleseismic  distances.  The  excitation  of  LONGSHOT  Rayleigh 
waves,  however,  was  an  order  of  magnitude  smaller  than  the  Rayleigh 
wave  excitation  for  an  equivalent  magnitude  earthquake  in  the  same 
region. 

5.  Love  waves  could  not  be  detected  for  the  LONGSHOT  explosion. 

6.  Teleseismic  surface  wave  signals  could  be  detected  for 
small  earthquakes  of  body  wave  magnitude  less  than  5  in  each  of  the 
three  different  source  regions  investigated.  The  actual  lower  limit 
for  each  region  is  still  uncertain. 
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Figure  1.  Schematic  Map  Showing  Source  Region  R  and  Trajectorie 
of  Surface  Waves  (of  a  Particular  Frequency)  from  Events  A  and  B 
to  an  Observation  Point  S.  Region  I  has  a  Higher  Phase  Velocity 
than  Region  II  in  this  Hypothetical  Case. 
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Figure  3A.  Synthetic  Signal  Test  Case 
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Figure  3C.  Synthetic  Signal  Teat  Caae 
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Figure  5A.  Effect  of  Eptcentral  uletance  Perturbation* 
on  Mat chad  Filter  Reaponae. 
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Figure  6.  Effect  of  Differences  in  Epicentral 
Distance  on  Correlation  Coefficient.  Dispersion 
Curve  Used  is  Shown  in  Figure  4. 


Figure  7.  Spectra  Used  to  Test  Amplitude  Effects  on  Matched 

Filter  Response. 
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Figure  8A.  Effect  of  Amplitude  Spectra  on  Matched  Filter 
Response .  The  Spectra  Ai  (f)  are  Shown  in  Figure  7.  The 
Diatance  ia  3000  km.  in  Each  Caae  and  the  Diaperalon  ia 
Shown  in  Figure  4. 
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Figure  9.  Change  In  Matched  Filter  Parameters  Due 
Shifts  in  Period  of  Spectral  Peak.  (t  «*  20  sec  ) 
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Figure  31.  Matched  Filter  Search  tor  LONGSHOT  (m 
Love  Wrvco  with  22  November  19&5  Earthquake  (m  -  5 
Love  Waves. 


Figure  12.  Observed  Values  of  e  with  Distance  for  LONGSHOT 
Relative  to  22  November  1965  Rayleigh  Waves.  Stations  are 
Along  &  Single  Azimuth  from  the  Source 


Figure  13.  Comparison  of  Rayleigh  Waves  from  a  Magnitude  4 
Event  (25  March  1965)  with  Those  from  a  Magnitude  5.8  Event 
(22  November  1965)  Near  Andreanof  Island. 
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Figure  J4A.  Search  for  Rayleigh  Wa-'es  From  a  Magnitude 
4.0  Even;-  at  Depth  48  Km.  Using  22  November  19r,5  Magni  ¬ 
tude  5.9  Event. 
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Figure  15A.  Matched  Filter  Comperieon  of  Rayleigh  waves  From 
a  Megnitude  5.0  Event  (29  January  1965)  with  Thoae  from  a  Ma7 
nit ude  5.3  Event  (20  April  1965)  ir.  Kamchatka. 
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Figure  16.  Observed  Relative  Excitation  (a)  of  Rayleigh  Waves 
Between  a  Magnitude  5.8  Event  (29  January  1965)  and  a  Magnitude 
5.3  Event  (20  April  1965)  Located  in  Kamchatka.  A  LONGSHOT  Pro 
file  is  Also  Shown. 
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Figure  18B.  Matched  Filter  Comparison  of  Rayleigh  Wavea 
From  a  Magnitud  5.8  Event  <29  January  1965)  With  Thoae 
From  a  Magnitude  5.0  Event  <14  February  1965)  in  Kamchatka 
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January  1965. 
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Figure  21B.  Search  for  Rayleigh  Wave*  From  a  Magnitude  4.1 
Event  (13  August  19M)  Ue:‘ng  Those  From  a  Magnitude  5.3 
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APPENDIX  I 


INDIRECT  METHOD  FOP  OBTAINING  SPECTRUM  OF  SMALL  EVEr^S 
From  the  analysis  in  the  text,  *c  'nave 

a(r>  n-i) 


The  frequency-domain  equivalent  of  ( I  —  1 )  is 

-  Co 

&(*)\  Y(w)TV)^  =  I  XMY^/e'^dv 

*  9  (1-2) 
where  Y(  lO  )  and  X(  (4. >  )  are  the  Fourier  transforms  of  y(t)  and 

X(  t)  respectively.  Without  losing  generality,  X  can  be  taken  as 
zero,  since  it  is  arbitrary  in  the  analysis. 

If  this  is  done,  we  have 


a  j  Y(uQY%)Ju>  ~  J  XMY 


(1-3) 


We  know  in  addition  that  in  propagating  over  a  distance 

A  from  the  spectral  amplitude  is  reduced  by  a  factor 

L 

r  *1*1  a 


where  2g  (  W  )  V  <«>/£<£  ,  Q  is  the  quality  factor  characteristic 

of  the  attenuation  due  to  anelasticity ,  u  is  radian  frequency,  and 
U (  U)  )  is  the  group  velocity.  We  sec  that  the  exponential  factor 


describes  the  amplitude  decay  due  to  departure  from  perfect  elasticity. 


and  the  term  jn  Luackctu  de&ei  i  ben  the  c-ft'ecl  of  gcomett  ical 

spreading . 

We  see  that  if  A  is  de terrr.i ned  at  scvcial  distances  alonu 


a  given  a;*-irnut5i  i  j  am  the  uoutcc»  we  ran  generalise 


(u), 


(1-3)  to 


2 


The  integrals  in  this  expression  are  now  in  the  form  of  ?  **tace 
transform  with  respect  to  A  .  Carrying  out  the  inversion  of 
both  sides  of  (1-6)  and  changing  the  order  of  integration  on 
the  left,  we  have 

.  J*ico  i 

c  [  S(ufc>**)\  a.(S)e 

J  C  ■'f-ito 


--  c  X(  sic, *,)Y* (*/*>**) 


(1-7) 


or 

c-‘J*Y*(y/c,  A.)  Z(i’U>')  £<*>'  5  C*X(5/c/Aa)y *(*/C,Ac) 

0  (1-8) 

That  is,  the  left-hand  side  is  the  convolution  of  Y  (s/c)  with 
A(  s) ,  the  inverse  transform  of  &(*)■  From  (1-8)  we  can 

find  X(s/e,  &o  )  directly  by  carrying  out  this  convolution  and 

* 

dividing  the  result  by  Y  (s/c  ,  A,  ).  Equivalently,  if  we  assume 
that  the  spectrum  of  x( t)  can  be  written  as  F(  u>  )Y(  ce  )  ,  then 
from  (1-8) , 

F  (  s/c ,  )  —  f  Y  (1-9) 

where  everything  on  the  right-hand  side  can  be  calculated  from 
observations  of  y(t)  and  A  (  ij  ) .  Then  the  spectrum  of  the  unknown 
event  at  distance  is  given  by 

«=  F(  Wj&t  )Y(  U»,  A0)  . 


3  - 


\ 

9 

i 


1W 


./a  >mi 


Results  for  a  particular  form  of  a.  (A) 

Suppose  that 

a,  '  \)  &  a  A 

where  G  and  q  nre  determined  from  a  least-squares  fit  to  the 
observations  of  CL  (  &  ) •  Then  using  this  approximation  in 
the  left-hand  side  of  (1-7),  we  get 

Qe'jVKsM  c ** 4u- 

,  rtc  MV(V/c,A<>)  if.-f-w'.M*'' 


Thus 


X(  0>,A*)  =  GY2[  &o]/YC&,  (1-11) 


Thus  we  can  calculate  the  spectrum  X(  A0),  since  we  can 
evaluate  the  right-hand  side  of  (1-11).  As  a  check  on  the 
result  (1-11),  we  see  that  if  there  is  no  change  in  Q,  with 
distance  (q=0) ,  then  X(  <M ,  be  )  =  GY(  ujf  ,3# )  =  a  (Aj)  Y (**>&*) 
implying  that  the  spectra  have  the  same  shape  -  or  equivalently, 
that  their  waveforms  are  identical.  If  there  were  no  attenuation 
with  distance  (c  =  0),  then  the  method  would  fail,  as  can  be  seen 
by  setting  c  =  0  in  (1-6) . 

In  practice,  estimates  of  g(  W  )  would  be  obtained  from  the 
observed  slope  of  log  [Kws&ty/Yfa,*)]  .  Then  the  best  linear 
approximation  (£U»)  to  g(  Lf  )  could  be  determined.  This  would  be 


1 


4 


the  preferable  approach,  since  the  total  propagation  losses  comprising 
both  intrinsic  absorption  and  scattering  would  be  accounted  for. 

In  summary,  then,  using  equations  (1-8)  or  (1-11)  we  have 
a  means  of  estimating  the  spectrum  of  the  unknown  event  in  terms 
of  the  spectrum  of  y(t). 


APPENDIX  II 


SURFACE  WAVE  SYNTHESIS 


We  have  written  a  program  which  computes  a  synthetic 
surface  wave  seismogram  from  assumed  group  velocity  and  source 
radiation  spectrum.  The  input  parameters  ares 
T  Period  in  seconds. 

C(T)  Phase  velooity  as  a  function  of  period. 

A(T)  Source  amplitude  spectrum. 

R  Epicentral  distance. 

Tx  Starting  time  (R/U^^). 

DT  Time  increment  for  synthesis  output. 

NP  Desired  number  of  output  points  in  synthetic 


seismogram. 


DF 

The  synthesis 


Frequency  increment, 
is  given  by: 


2 


Cos 


(it-.; 


Ah'  Mfl 


}  ] 

II 


U>i  ft- R/C(w/)J 

Avlt-R/Vc»0] 


where 


(group  velocity) 


r  u>»  fdC(u>A 

Vl*i)  '  C(i»*) I  C(ml)  \  Au>  )u 

This  is  a  generalization  of  the  synthesis  method  devised  by 
Aki  (Reference  2). 

For  a  derivation  of  Equation  (II-l),  see  Reference  3,  pp.  139. 
We  convert  input  periods  to  frequency,  and  fit  polymonials  to  the 
input  phase  velocity  and  amplitude  values,  both  to  facilitate  the 
calculation  of  U(  W  )  and  dAfdv  ,  and  to  allow  us  to  evaluate  the 
integrand  in  Equation  (II-l)  at  equal  increments  of  frequency. 
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